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The densest local packing (DLP) problem in d-dimensional Euclidean space R*^ 
involves the placement of N nonoverlapping spheres of unit diameter near an 
additional fixed unit-diamctcr sphere such that the greatest distance from the 
center of the fixed sphere to the centers of any of the N surrounding spheres is 
minimized. Solutions to the DLP problem are relevant to the realizability of pair 
correlation functions for sphere packings and might prove useful in improving 
upon the best known upper bounds on the maximum packing fraction of sphere 
packings in dimensions greater than three. The optimal spherical code problem 
in M"' involves the placement of the centers of N nonoverlapping spheres of unit 
diameter onto the surface of a sphere of radius R such that R is minimized. It 
is proved that in any dimension, all solutions between unity and the golden ratio 
T to the optimal spherical code problem for N spheres are also solutions to the 
corresponding DLP problem. It follows that for any packing of nonoverlapping 
spheres of unit diameter, a spherical region of radius less than or equal to r 
centered on an arbitrary sphere center cannot enclose a number of sphere centers 
greater than one more than the number that can be placed on the region's surface. 



I. INTRODUCTION 



The densest local packing (DLP) problem in seeks an arrangement of N 
spheres of unit diameter near (local to) an additional fixed central sphere such 
that the greatest radius R between the centers of the surrounding spheres and 
the center of the central sphere is minimized. For an optimal configuration of 

spheres, i.e., a configuration for which R is minimized, we call the minimized 
greatest radius R^^^i^N). The 'Z' in the notation R^^^{N) serves to distinguish 
from i?^-„(A^), where -Rmm(^) optimal spherical code (OSC) problem is 

the radius of the minimal radius sphere onto the surface of which can be placed 
the centers of nonoverlapping spheres of unit diameter. For = 15, c? = 2, Fig. 
[l] depicts a conjectured optimal configuration for the DLP problem with minimal 
radius i?^j„(15) = 1.873123 . . . alongside an optimal spherical code configuration 
with minimal distance i?;^j„(15) = 2.404867. . . . 




FIG. 1. Left: a conjectured optimal DLP configuration for N = 15, d = 2, 
i?^j„(15) = 1.873123 .... Right: a spherical code optimal configuration for N = 15, 
d = 2, i?^i„(15) = l/(2sin(7r/15)) = 2.404867.... 

The kissing number problem in M"^ seeks the maximum number of nonover- 
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lapping spheres that may simultaneously be in contact with a (additional) 
spheref^ it is a special case of the DLP problem in that is equal to the 
greatest for which Rfaini^) ~ ^- '^^^^ DLP problem can also be said to en- 
compass the sphere packing problem in that in the limit as — t- oo, optimal 
sphere packings and optimal DLP packings are equivalent. 

The maximum possible with R = R^^^{N) for an optimal DLP configura- 
tion of spheres in M*^ is the maximum of the function Z{ri,R). The function 
Z(rj, R) is defined for packings of nonoverlapping spheres of unit diameter as 
the number of sphere centers that are within distance R from a sphere center at 
position Fj, with i an index over all centers and where the value of Z{ri, R) does 
not count the sphere center at Fj. For a statistically homogeneous packing, the 
maximum at fixed R of Z{ri, R) is an upper bound on the maximum of the func- 
tion Z{R), where Z{R) is the expected number of sphere centers within distance 
R from any given sphere center, or equivalently the average of Z{ri, R) over all i. 
For a packing that is also statistically isotropic, Z{R) can be related to the pair 
correlation function g2{r), a function proportional to the probability density of 
finding a separation r between any two points and normalized such that it takes 
the value of unity when no spatial correlations are present, by 



Jo 

where p is the constant number density of points and Si(r) is the surface area of 
a sphere of radius r in d dimensions. 



The optimal spherical code (OSC) and DLP problems are similar. A spherical 
code is defined for parameters {d, N, t) as a set of A^ vectors from the origin 
to points on 

gd-i ^ ^^^Yi that the inner product between any two distinct 






(2) 



3 



vectors is less than or equal to t. The OSC problem is to minimize t given A^, 
or to maximize N given t. There have been a number of investigations into the 
optimality and uniqueness of specific spherical codes (for example, see Ericson 
and ZinovieA^ and Cohn and Kumail^, and into providing bounds on given t 
and 

A spherical code may be represented by a packing of nonoverlapping spheres 
of unit diameter with centers distributed on the surface of a sphere of radius R. In 
this representation, the OSC problem for a given A^ requires finding the minimum 
R, R^^^{N), such that no two spheres overlap, i.e., such that the distance between 
the centers of any two spheres is greater than or equal to unity. 

The OSC problem formulated in terms of nonoverlapping spheres and the DLP 
problem differ for all A^ where i?,^j„(A^) > 1 only in that the former restricts the 
placement of sphere centers to a subset of the space allowed in the latter. From 
this observation, it is clear that when there exists a configuration of spheres that 
is a solution to the DLP problem with minimal radius Rmmi-^) ^^^^ ^^^^ ^ 
spherical code, it is additionally a solution to the corresponding OSC problem, 
with RL^{N) = Ri,^{N). 

II. THE DENSEST LOCAL PACKING PROBLEM AND 
REALIZABILITY 

Only functions obeying certain necessary conditions known as realizability 
conditions can be correlation functions of point processes in M'^.HHSl Two realiz- 
ability conditions on the pair correlation function g2{r) are the nonnegativity of 
g2{r) and its corresponding structure factor S{k) at all points r and fc.'^ These two 
conditions appear to be strong conditions for the realizability of sphere packings 
(point processes in which the minimum pair separation distance is unity), espe- 
cially as the space dimension increases.!^ They have been employed, among other 
uses, to provide conjectures for a lower bound on the maximum packing fraction 
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of an infinite sphere packing in any dimension, and to demonstrate the feasibihty 
in three dimensions of a sequence of disordered packings whose disorder vanishes 
as density approaches the maximum possibleP 

Cohn and Elkied^ employ analogs of these two conditions, in conjunction with 
a linear programming technique, to find the best known bounds on the packing 
fraction of infinite sphere packings in (at least) dimensions four through 36. In 
the conclusions of a previous work,- we discuss how a third realizability condition, 
found by solving the DLP problem for a packing of 13 spheres in three dimensions, 
can improve upon the three-dimensional bound found in Ref. 9. 

The technique employed in Ref. 7 to find conjectured lower bounds has been 
shown to be the dual of the primal infinite-dimensional linear program employed 
in Ref. 9, and Cohn and Kumai'^have since shown that there is no duality gap 
between the two programs. This means that when the best g2{r) test functions 
are employed, the upper and conjectured lower bounds will coincide. Cohn and 
Elkies in Ref. 9 were able to find a test function that yields the best upper bound 
on the maximal packing fraction in three dimensions, a packing fraction of 0.778, 
which is well above the true maximum. This means that there is a test function 
for the lower bound formulation that will deliver the same packing fraction of 
0.778, which is clearly not realizable. 

A putative improvement on the upper bound in M.^ was obtained by employing 
an estimate for i?^j„(13) in the DLP problem in M^.'m Requiring that Z{R) < 
12 up to some small positive a beyond contact, with R = 1 + a the estimate 
for _R^^„(13), reduces the d = 3 bound in Ref. 9. For example, estimating 
a = 0.0^^ reduces the bound from 0.778 to 0.771. This result strongly suggests 
that DLP solutions introduce more information than is contained in the pair 
correlation function alone, in that there is at least one test g2{r) that obeys the 
two nonnegativity conditions but violates the bound Z{1 + a) < 12. 

Further solutions to the DLP problem provide additional realizability con- 
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ditions that might be employed to improve upon the upper bounds on infinite 
sphere packings in dimensions greater than three. For a statistically homogeneous 
and isotropic packing of spheres, these additional conditions may be written as 

Z{R) < Z^a.{R), (3) 

where the function Z^axiR) is defined in as the maximum number of sphere 
centers that fit within distance R from a central sphere center .1^ It is clear that 
Zmax{R) in is completely defined by the solutions to the densest local packing 
problem at all A^. 

In the following section, we show that any configuration of d-dimensional 
spheres near a (additional) sphere fixed at the origin, with the greatest of the N 
distances from the origin to the N sphere centers equal to i? < r = (1 + -\/5)/2 ^ 
1.618034 the golden ratio, may be transformed to a spherical code in the sense 
of nonoverlapping spheres, also of radius R. As this statement is applicable 
to any configuration of A^ spheres that is a solution to the DLP problem with 
R < T, it follows that any optimal spherical code with radius i?,^j„(A^) < r is also 
an optimal configuration for the corresponding DLP problem, with i?^j„(A^) = 

RLniN). 

III. TRANSLATING UNIT-DIAMETER SPHERES TO THE 
SURFACE AT RADIUS R<t 

The key idea behind the proof of the above statement involves translating 
sphere centers radially outward to a spherical surface of radius R. The idea of 
radially translating points to a spherical surface has been employed by MelisserP^l 
to aid a proof of the optimality of certain packings of 11 congruent nonoverlapping 
circles in a circle and more recently by Cohn and Kumai^ to rescale vectors in 
to terminate on 5*^^. However, prior to this work, the maximum radius from 
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the center of a fixed nonovcrlapping sphere to which the centers of surrounding 
spheres can be translated without resulting overlap was not known. 

Specifically, for any number of nonovcrlapping spheres of unit diameter ini- 
tially situated such that their centers are contained in a spherical shell of radial 
span [1, R] with 1 < i? < r, all sphere centers at a distance less than R from 
the center of the shell may be translated radially outward to distance R without 
any resulting overlap between spheres. This statement more generally applies 
(via a simple scaling argument) to congruent nonoverlapping spheres of arbitrary 
diameter D that are contained within a spherical shell of radial span [D,R], 
D <R<tD. 

Define An{R) in R** as the set of all packings of any number N of nonoverlap- 
ping spheres of unit diameter with centers situated in a spherical shell of radial 
span [1, R], with R the greatest of the distances from the center of the shell (the 
origin) to the N sphere centers. An element of the set Ais[{R) represents any 
arrangement of N spheres with greatest distance R situated near an additional 
nonoverlapping sphere fixed at the origin. 

Theorem 1. Consider any single element of An{R) in R'^. For R < r, all N 

spheres may be translated radially outward such that their centers are at distance 
R from the origin and still remain an element of Ai\f{R). For R > t , d > 1, there 
exist elements of Aiy{R) such that an outward radial translation of a given sphere 
center to distance R will yield overlap between at least two of the N spheres. 

Proof. The proof proceeds from the law of cosines in the method of the proof of 
Lemma 4.1 in Ref. 15. For any two of the N spheres with centers situated at 
distances b, c from the origin and separated by distance a, the cosine of the angle 
formed between the two centers at the origin, taken such that < 9 < n, is 



For nonoverlapping spheres of unit diameter, a > 1, and 

^"^^^ 26c ' 
where the equahty holds when the two spheres are in contact. 

Over the range h> 1, c > 1, the function cos^^ in Eq. ([s]) is convex individually 
in both h and c. This implies that cos^ must be at a maximum at one of the 
corners of the square l<h<R, l<c<R. If R > r, the point (1, R) (or 
equivalently (-R, 1)) yields the maximum, whereas for R < r, the point {R,R) 
yields the maximum, with (1, R) and {R, R) both yielding the maximum at R = t. 

It follows directly that for R < t, the minimum possible angle at the origin 
between any two of the centers of spheres that are an element of Ajv(i?) is the 
angle present when two of the centers are placed at distance R from the origin 
and distance unity from one another. An outward radial translation of one or 
both of any pair of centers to distance R will therefore yield no overlap between 
the two spheres, as the angle between the centers must be greater than or equal to 
the angle present when two spheres are in contact with one another with centers 
at distance R from the origin. As this holds for any pair of the A^ sphere centers, 
all centers at a distance less than R, R < t, may be translated radially outward 
to distance R without any resulting overlap. 

For R > T, d > 1, overlap between two spheres is possible after an outward 
radial translation. For example, when two spheres are initially in contact with 
centers at distance unity from each other and at distances 1 and R> t from the 
origin, the angle formed at the origin between centers is smaller than the angle 
present when the spheres are in contact with centers both at distance R. A radial 
translation outward of the sphere center at distance 1 to distance R would thus 
yield overlap. This concludes the proof of Theorem 1. □ 
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IV. RESULTS AND DISCUSSION 



Theorem 1 applies to any configuration of nonoverlapping splieres tliat 
are an element of Afyf{R). In particular, in an optimal DLP configuration 
for splieres in M*^ with R^^^{N) < r, which is by definition an element of 
An{R^^^{N)), any of the spheres with centers not at distance R^^^^N) from 
the origin may be translated radially outward to distance Rf^^^{N) without any 
overlap between spheres. The resulting configuration is both a solution to the 
DLP problem and, in the sense of nonoverlapping spheres, to the corresponding 
OSC problem, with R^ini^) — ^mini-^)- Theorem 1 therefore implies that while 
for 1 < R^^^{N) < T there may be solutions to the densest local packing problem 
that are not spherical codes, for 1 < R^^.^{N) < r there are no solutions to the 
optimal spherical code problem that are not also solutions to the corresponding 
densest local packing problem. 

The kissing numbers in M"^ are only known rigorously for d = 1 . . 4, d = 8 
and d = 24^foT = 1, 2, 3 and 4, they are 2, 6, 12 and 24,IIS1 respectively. For 
N < Kd, the solution to the DLP problem is simply R^^^{N) = 1 by necessity 
as the nonoverlapping sphere of unit diameter at the origin is fixed. For such 
that Kd < N < N^, where we define iVJ in M'^ as the greatest integer N such 
that R^i^{N) < r, the optimal spherical codes are solutions to the corresponding 
DLP problems with Rmini-^) — -^mm(^)- The questions concerning the values 
of in each dimension and how grows with d naturally emerge. 

In one dimension, the answer to the first question is trivial, with A^^ = 2. 
In two dimensions, optimal spherical codes can be found analytically via simple 
trigonometry, with R^^^{N) = r ior N = 10, or A'J = 10. Strong conjectured 
solutions for R^^^{N) that serve (at least) as upper bounds to the OSC problem 
are well-known in low dimensions greater than two for small A^.f^^For c? = 3, these 
yield the conjecture = 33 with R^-^{33) ~ 1.607051. For d = 4, a. unique 
optimal spherical code is known such that i?^j„(120) = r, giving the result that 
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Nl = 120P 

The question of precisely how grows with d is still open, and is more 
complicated; however, bounds may be established via known bounds on (given 
d and t) for optimal spherical codes, such as with those given in chapter two of 
Ref. 1. The lower bound (due to Wyner^-) on N[d, 0) for a spherical code of 
minimum angle = cos^^(t) in dimension d is 

iV(rf,0)>^47TV' (6) 
sm (0) 

giving for R^^^{N) = r, t = r/2 and > I.TOIS'^. This may be compared to 
the lower bound on the kissing number obtained from ([6|, > 1.1547'^. The 
upper bound due to RankirP^ is 

Nid,<P)<(-nd'cosi<P)j (v^sin(0/2)) , (7) 

giving, for R^.^W = r, < 1.12736/3/^x2.2883'^. This may be compared to the 
Kabatiansky-LevenshteirpH upper bound on the kissing number, Kd < 1.3205'^. 
Comparing the upper bound on the kissing number and the lower bound on N^, 
it is clear that grows exponentially faster than K^. 
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